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A supersymmetric extension of the Standard Model, witli maximal sn-bR mixing and a new 
source of CP violation, contain all the necessary ingredients to account for a possible anomaly in 
the measured CP asymmetry in i? — > <j)Ks decay. In the same framework we study the decay 
B J/tpK*, paying particular attention to observables that can be extracted by performing a 
time dependent angular analysis, and become nonzero because of new physics effects. 
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I. INTRODUCTION 

Recent results on B decay to charmonium modes, such 
as i? — *■ J/ip Ks{K*), are in good agreement with the 
Standard Model (SM). However, for the time-dependent 
CP asymmetry in S — > (pKs mode, also if new results 
seem to show more agreement with the SM, S^Ks = 
0.50 ± 0.25 (BaBar) and S^Ks ^ 0.06 ± 0.33 (Belle) 
0, than in the past, (before ICHEP04, the world average 
S^Ks = 0.02 ± 0.29 [|), New Physics (NP) effects could 
still play some role. A deviation from the SM expectation 
of S^Ks — S-ifiKs would call for large s — b mixing, the 
existence of a new source oi CP violation, and perhaps 
right-handed dynamics 0] ■ It would be important to find 
confirming evidence in the future. 

In anticipation of a future "super B factory" that 
would allow precision measurements, we study CP viola- 
tion in the vector- vector J/tpK* mode by taking 
into account deviations in ScfiKs- We do not expect the 
NP effects to show up in the branching fraction since, in 
contrast to B ^ <t>Ks, B^ J/ipK* decay is tree dom- 
inant. With special attention to observables related to 
the time dependent CP asymmetry 0,1^13) we focus on 
manifestations of NP effects. In the phenomenological 
analysis presented in this paper, we pay particular atten- 
tion to those obserables which are expected to vanish in 
the SM 0. 



the amplitude is given by changing the sign of the weak 
phase (j). One then obtains the direct CP asymmetry 
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which does not vanish for ((> ^ and if the strong phase 
difference is also not zero. However, since the b quark is 
rather heavy, the strong phases are expected to be small 
and the effect of NP in the direct CP asymmetry could 
be washed out. It is therefore important that one can 
still seek for NP effects by performing a time dependent 
analysis and comparing (i) J/tpK* with. J /tpKs- In 
fact, more information is contained in the time dependent 
angular analysis of vector- vector decays such as B^{t) 
J/iPK* or (t)K* 

For a B ^ VV decay, the final state can be decom- 
posed into three helicity amphtudes {Aq, A|| , A^}. Ao 
corresponds to both the vector mesons being polarized 
along their direction of motion, while A|| and A± cor- 
respond to both polarization states being transverse to 
their directions of motion but parallel and orthogonal to 
each other, respectively 8j. If in particular we consider 
the decay B^{t) J /ipK* , analogous to Eq. we have, 
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Ax{B^J/^K*) = axe' 



bxe^"t'e'^^ 
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The J/^pK^*'> decay is dominated by the tree 

level b — > CCS process, while the CP phase in the corre- 
sponding penguin amplitude is highly suppressed. If NP 
contributions are present, it could manifest itself as di- 
rect CP violation effects. This can be illustrated by the 
full amplitude for the decay B f, 
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be'^e'^ 
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where the weak phases are assumed to be zero for the first 
amplitude and (p for the second, and (5°'^ are the respec- 
tively strong phases. For the CP conjugate decay B f 



where ax and bx are the SM and NP amplitudes and S"^'^ 
their respective strong phases, for each helicity compo- 
nent. The full decay amplitude becomes 

A{B ^ J/^K*) = Aogo + A\^g^\ + lA^g^, (6) 
AiB ^ J/^K*) = Aogo + - lA^g^, (7) 

with gx the coefficients of the helicity amplitudes in the 
linear polarization basis |^. If one considers the case 
where K* and K* are detected through their decay to 
KsTT^ so that both i?° and 5° decay to a common final 
state, the time dependent decay rates can be written as. 
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By performing an angular analysis and time dependent 
study of the decays i?^ — + J/ipK* and — > J/^K*, one 
can measure the observables A\„ , Eact and g\a- • These 
observables can be expressed in terms of the normalized 
helicity amplitudes ^o, ^ii and Aj_: 



where i = {0, ||}, q/p — exp {—2i(f>^i-x) with 0,nix the 
weak phase in B^ — B^^ mixing. From Eqs. Q and © 
one can obtain the same observables in terms of a\, b\, 
(j>, Sx = 5^- Si and = (5^^ - 7]. In particular, A^, 
can be expressed as, 
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-Im ( -i^l^j 



lm[^iAlAo + A*A^\) 
Im ( ^A*A, , , 



A±i = 2[aj_6i cos(Ai — — a^foj^ cos(Ai+(5^)] sin(/) . (10) 

The observable Aj^^ is special, as made clear by Eq. pU|l . 
because it remains nonzero in the presence of NP effects 
{(j) ^ 0), even if the strong phase differences vanish. In 
contrast, direct CP asymmetries E^a are washed out if 
the strong phase differences vanish 



(9) 



III. B^J/ipK* 



We can now proceed towards NP effects in i? ^ 
J/ip K* . We start by writing the decay amplitudes for 
B — > J/ijj K* using the factorization approximation, but 
keeping the color octet contribution llf 



Ax{B~^J/^K*) 



i^VcbV:j^m^e';{X){af{K* \ sj,{l - 75)6 I B) 
as mb , 
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where A = A 



J/V 



Ak* = 0, ±1 denote the hclicities 



of the final state vector particles J/i/j and K* in the 5° 
rest frame [lll |. The dominant contribution in Eq. Hll() 
is given by the tree level term proportional to 
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while the color dipole moment terms, with the C12 oper- 
ator coming dominantly from SM and the c'12 operator 
due exclusively to NP, give smaller corrections. In the 
expression for a|*^ we have neglected the strong and elec- 
troweak penguin contributions. The quantity C — j^+Cs 
in Eq. H12I) takes the value 1/Nc = 1/3 in the naive 
factorization while deviations from 1/A'^c due to non- 
factorizable contributions to the hadronic matrix ele- 
ments are measured by the parameters and |lO| . 



The effective Wilson coefficients c^j^jmp) and C2^(mb) for 
a 6 ^ s transition are defined in Ref . [ij ■ 

The way we proceed to determine the parameter C fol- 
lows Ref. ^3 . We fit the branching ratios for the decays 
B J/tPKs{K*) and B ■ipi2S) K* to extract af . 
We checked explicitly that the NP effect of Ref. Q does 
not make significant impact on the decay rates. However, 
the extraction of a'^ depends on the specific model one 
uses for the hadronic form factors ^13.] . In this work we 
use the form factors at zero momentum transfer for the 
B V transitions obtained in the light-cone sum rule 
(LCSR) analysis The form factor dependence is 
parametrized by 
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the squark mass scale, and 
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FIG. 1: c'i2 as a function of a for mi = 200 GeV and m = 
2 TeV. Re (c'12) and Im (c'12) are plotted with solid, dash and 
solid-dash lines for mg = 300, 500, 800 GeV. 



where the values of the parameters a and h are given in 
Ref IT3|. Our extracted value for is 0.20, which using 
Eq. (fT^ and the effective Wilson coefficients of Table 1 of 
Ref. (12i gives C, = 0.48, or the effective number of colors 
1/C = 2.1. Knowing C, one derives = 0.15 and for the 
value of we assume = . 



IV. B J/itK* AND NP EFFECTS 

As previously mentioned, in this work we focus on the 
decay B'^{t) J li\)K* in the context of a supersymmet- 
ric model with maximal sr — 6/? mixing. An approximate 
Abclian flavor symmetry |l5l | can be introduced to justify 
such a large mixing. In this model the right-right mass 
matrix m"^^^^ for the strange-beauty squark sector takes 
the form 



M 



2{sb) 
RR 



'"■22 '"'23^ 



33 



R 



R\ 



(14) 



The phase a is the NP weak phase which will affect the 
observables one can extract from an angular analysis. Be- 



cause of the almost democratic structure of M 



2(sb) 
RR ' 



one of 



the two strange-beauty squarks, sbi, can be rather light, 
even for rh ~ O(TeV). A light strange-beauty squark, 
together with a light gluino, can make S'^^/f^ negative for 
cr < 7r/2 ||. 

The main new contribution to -B°(t) — > J/tpK* is 
given by the color dipole moment amplitude through 
gluino and sb squark exchange in the loop. The an- 
alytic expressions for the Wilson coefficient c'12 for 
the color dipole moment operator q., / (871^)771^ s„ g^'" ( 1 ± 
gz)\^p/2 bpG^^ can be found in Ref. 

C12 remains basically the same as 
we plot the real and imaginary parts of the Wilson coef- 
ficient for three different values of gluino mass mg = 
300, 500, 800 GeV The eigenvalues of Eq. (0 are 
taken (with some level of tuning) as mf = (200 GeV)^ 



,/i87T 



'12 — 



The coefficient 
0.15. In Fig. 1 



and mi = 2rh 



^ 2 



with = (2 TeV)" 



the 



Using the parametrization for 
{K* I S7^(l -75)6 I B) and {K* \ sia, 
given m Ref. [ill the amplitudes Ax{B 
be written as. 



matrix elements 
^,7m(1±75)^> I B) 
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C8(ci2 + €[2) 9+{ml)e'"''''^ pK'aPBp + i^^£.s{ci2 - ci2) 



2tt ml 
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(16) 



The general covariant form for Ao,±i is given by, 
^o,±i =e;^(0,±l)ejf.,(0,±l) as'' 



■ PbPb + 



m^mx 



-e^^'^^PK'cPBfS 



(17) 



where a, b and c are three invariant amplitudes. The corresponding amplitudes Aq ±i are obtained by taking the 
conjugate of the invariant amplitudes a, b, c and switching the sign of the term c^'^"^ PK-aPBi3 in Eq. p7|) . By 
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comparing Eq. H16() with Eq. H17() one can extract the three invariant amphtudes 

as m 
27r m 



off . Q^s f. I I \ ( t 2\ j,/2n2\ 



6 = —m^mK' 



tL 



V{ml) 



rriB + niK' 



off OLsms ^ , I , ,2 . 



(18) 
(19) 
(20) 



In the expressions of Eqs. p8|l - H2()|l we have omit- 
ted the common factor \/2Gi? Vc^V^* /^m^. Note that 
each invariant amphtude contains a SM contribution 
which is dominated by the tree level term proportional 
to a'^ plus the color dipole moment term proportional 
to C12, and a NP contribution proportional to c\2- One 
can consequently write the three invariant amplitudes 
as the sum of a SM and a NP contribution: a, 6, c = 
(a, 6, c)SM + (a, 6, c)^?. 

We now rewrite the helicity amplitudes Ao.±i in terms 
of the invariant amplitudes and the kinematic factor x = 
Pip ■ PK' I {jn^rriK') 0|. Writing separately for the SM 
and NP contributions, we obtain, 



SM,NP 
±1 

SM,NP 




^SM,NP _|_ ^SM,NP^ 



.^SM,NP^ 



X 

^SM,NP 



2 _ 1 
2 



(x^-1), (21) 



with ^o,±i = -^oXi +^a±i- To evaluate the observables 



of Eq. ^ one can transform ^q'^';^'^ in the correspond- 
ing linear polarization amplitudes ^q'||''^^ using the re- 
lations: ^11, _L = (^+1 ± ^-i)/\/2, ^0 being the same in 
both basis. 

The invariant amplitudes in the linear polarization ba- 
sis ao,||,-L and foo,||,-L can subsequently be expressed in 
terms of (a, 5, c)'(sm,np)^ 









= -aNPx-6NP(^2 


ail 


= V2aS^, 




= V2aN^ 


a±_ 


= c^^^2{x^ - 1) , 




= c^V2(:.2-l), 



1), 
1), 



(22) 

where (f> = arg(c']^2) ^ tt. 

In Fig. 2 we plot the observables K±i versus the 
NP weak phase a = arg(c']^2) for three different val- 
ues of the gluino mass nig = 300, 500, 800 GeV, with 
squark masses m\ = (200 GeV)^ and = 2to^ — 
with m = 2 TeV. We note that the effects of NP can 
be at most a few percent, and tend to disappear as 
the gluino becomes heavier |22j . This can be under- 
stood from Eq. (|10|l together with the expressions for 



contributions to A;^;^ are proportional to a\ (tree level 
dominated) hence are of 0(1), while for K±i one has 
0{b\/a\) ~ 0(0.01) at most, being suppressed by the 
ratio pNP^jiSM -^j^jj p g^j^j indicating respectively the 
penguin and tree terms. 

We see that, for the particular model considered in this 
work, to observe NP effects by performing an angular 
analysis for the decay -Bj](t) — > J/4'K*, one needs to 
be able to extract A±i with a precision of at least _a few 
per cent. On the other hand, as stressed in Rcf. 7], no 
tagging or time dependent measurements are needed to 
measure A±i since it appears with the same sign in both 
rates for B^{t) and B°i{t) (see Eq. ©). 

Following from the above considerations, it is evident 
that the decay B'^{t) —* (pK* becomes really interesting. 
This decay, contrary to B'^ (t) ^ J/ip K* , is not tree level 
dominated. Rather, it is of pure penguin type, and in the 
model considered a\ and 6a are of the same order. This 
implies now that the observables A±i not only differ from 
zero if there are NP effects, but they are expected to be 
of 0(1), being the ratio pNP/pSM _ q^-^-^^ Obviously for 
the decay B'^ (t) (pK* the hadronic uncertainties play a 
more important role than for the decay B'^ {t) J/ip K* , 
plaguing the theoretical prediction for Aj^^. Furthermore, 
the strength of the transverse components are not yet 
understood. 

For -Bj](t) — > J/ip K* decay, full angular analysis by the 
CLEO Collaboration shows that the P wave com- 
ponent is smaU, | P p = | p = 0.I6 ± 0.08 ± 0.04, 
while the longitudinal component is around 50%, Tl/T = 
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FIG. 2: Axi as a function of a for mi = 200 GeV and m = 
2 TeV. Axi are plotted in solid, dash and solid-dash lines for 
m-g = 300, 500, 800 GeV. 
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FIG. 3: I P p as a function of a for rrii = 200 GeV, m 
2 TeV and = 300 GeV compared with experiment. 



l-i 




sin 201 (S J/^Ks) in Fig. 6 with sin205"[f'^^ one 
can observe that deviations from sin 2(f)™^^ have opposite 
signs. ;25i This divergent behavior could in principle be 
easier to be observed than the single deviations. 



V. CONCLUSION 

A supersymmctric extension of the SM with a hght 
right-handed "strange-beauty" squark, a light gluino and 
a new CP phase, seems to contain all the necessary ingre- 
dients to explain the recent CP anomaly in Bd — > (j)Ks. 
In the same framework we have calculated possible NP 
effects to observables that can be extracted by the time 
dependent angular analysis of — > J/^p K* . An impor- 
tant role is played by the quantities Aj^^ with i — {0, ||}, 
which can be nonzero in the presence of NP even for very 
small strong phase differences. Our results show that de- 
viations from zero can be at most of the order of a few 
percent, since it is suppressed by the ratios of the NP 
penguin amplitude to the SM tree amplitude. This obvi- 
ously suggests that for decays which are pure penguins, 
like Bd — > 4>K* , deviations from zero for the observables 
K^i are expected to be of order one. 



FIG. 4: Vl/T as a function of a for rfii = 200 GeV, 
2 TeV and nig = 300 GeV compared with experiment. 



I Ao I 0.52 ±0.07 ±0.04. Recent measurements for the 
longitudinal and transverse amplitudes have been also re- 
ported by both BaBar and Belle |23| Collaborations 
with the respective values \A^\'^ = 0.16 ± 0.03 ± O.OI, 
I Ao 1^ = 0.60±0.03±0.02and | P = 0.19±0.02±0.03, 
I Ao 1^ = 0.62 ± 0.02 ± 0.03. In Figs. 3 and 4 we plot re- 
spectively I P 1^ and Fi/F as a function of the NP weak 
phase a. It seems clear that both statistical and system- 
atic errors need to be reduced by an order of magnitude 
to discriminate a nonzero value for Aj^^ as predicted by 
the model considered. We expect that the error on the 
extracted value of A^i will be of the same order as the 
one on I P 1^ or F^/F. 

We conclude this section by presenting the results for 
s\rL2(t)'^^'^ = Im(^|4l)' ^here A = 0, ||,± |2|. If NP 

effects are present one will not be able to extract (f)™™ 
[2^ (the phase of B^ — B^ mixing which in general can 
be affected by NP) and the measured value of which 
will depend on the helicity of the final state, will dif- 
fer from the real value of 0""'^ In Fig- 5 '^e plot 
sin 2(j)™^'^^ as a function of the NP weak phase a. As for 
the quantities A^^, the effects of NP on sm2(j)'^^^^ can be 
at most of few per cent. Deviations from sin 20™"' reach 
their largest value at cr = 7r/2 and are bigger for the 
transverse components, {A =||,±}. On the other hand 
deviations on sin20"Q'''*^, even if smaller, could be eas- 
ier to detect because of the higher number of longitudi- 
nally polarized final states. On top of that by comparing 
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FIG. 5: sin 205";; 



as a function of a for mi = 200 GeV and 
m = 2 TeV. sin 205°;f^= 



dash lines for mg 



are plotted in solid, dash and solid- 
300, 500, 800 GeV. The black solid line 



corresponds to sin 20]°"' = 0.733 
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FIG. 6: sin 201 (B — > J/tpKs) as a function of a for mi = 
200 GeV and m = 2 TeV. sin20i(B J/^jKs) is plotted in 
solid, dash and solid-dash lines for m^ = 300, 500, 800 GeV. 
The black solid line corresponds to sin 20™'" = 0.733 
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The quantities sin 2(j)™^'^^ can also differ from the real 
value sin 205""^ because of NP effects. In particular 
s'm2(j)flf'^ and sin20i(_B J/tpKs) have opposite de- 
viations from sin 205""^, and by comparing the two be- 
haviors, NP effects could be easier to be observed than 
by looking at single deviations. As for Aj^^, wc found 
that deviations are of the order of few percent. 

In conclusion, New Physics effects to Bd — > J/ip K* 
from the model considered in this work are found to be 
too small to be observed at the current B factories. But 
because of the small NP effects, B^(t) — > J /^p K* remains 



a good mode for measuring sin 20™"^. 
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